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The Lyapunov exponent spectrum and covariant Lyapunov vectors are studied for a quasi-one- 
dimensional system of hard disks as a function of density and system size. We characterize the 
system using the angle distributions between covariant vectors and the localization properties of 
both Gram-Schmidt and covariant vectors. At low density there is a kinetic regime that has simple 
scaling properties for the Lyapunov exponents and the average localization for part of the spectrum. 
This regime shows strong localization in a proportion of the first Gram-Schmidt and covariant 
vectors and this can be understood as highly localized configurations dominating the vector. The 
distribution of angles between neighbouring covariant vectors has characteristic shapes depending 
upon the difference in vector number, which vary over the continuous region of the spectrum. At 
dense gas or liquid like densities the behaviour of the covariant vectors are quite different. The 
possibility of tangencies between different components of the unstable manifold and between the 
stable and unstable manifolds is explored but it appears that exact tangencies do not occur for a 
generic chaotic trajectory. 
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I. INTRODUCTION 

The study of the Lyapunov vectors of particle systems 
generated by the Benettin scheme [1, 2 - the so-called 
Gram-Schmidt (GS) vectors - is now well advanced for 
quasi-one-dimensional systems [5 and two dimensional 
systems [3 . The existence of a stepwise structure in the 
smallest non zero Lyapunov exponents is an extensive 
feature of the Lyapunov spectrum (the full set of Lya- 
punov exponents) which has been well studied [SHSl [12] . 
For each step in the Lyapunov exponent spectrum the 
associated Lyapunov vector has a delocalised wavelike 
structure - referred to as a Lyapunov mode - which is of 
a particular type; either transverse, longitudinal or mo- 
mentum proportional [6^. These modes are connected 
to the slowest macroscopic fluid properties (the dynam- 
ically conserved quantities) or equivalently the trans- 
port of mass, momentum or energy. Transport of con- 
served quantities is always the slowest processes in par- 
ticle systems thus they relate to the smallest instabilities 
or smallest Lyapunov exponents. 

The statistical mechanics of chaotic many particle sys- 
tems is an illustrative example of a probabilistic treat- 
ment of the global behavior of deterministic microscopic 
dynamics. A great deal of effort has been devoted to 
finding links between macroscopic fluid quantities, such 
as transport coefficients, and chaotic properties of mi- 
croscopic systems such as the Lyapunov exponents [33]- 
[35] , There have been some significant successes such as 
the conjugate pairing rule for the Lyapunov spectrum in 
some thermostated systems jTOl [T3j [15] and the fluctu- 
ation theorem [Ml [16] yet there have also been illustra- 
tions of non-chaotic systems with well defined transport 
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properties. 

Microscopic particle systems that are models of real 
fluids can also be described from a macroscopic view- 
point. The description begins with conservation equa- 
tions and leads to the Navier- Stokes equations of fluid 
dynamics. The time evolution of the macroscopic quanti- 
ties in the Navier-Stokes equations can also be studied as 
a chaotic dynamical system but their phase space behav- 
ior is fundamentally different to the phase space behav- 
ior of the original particle system. Nonlinear, dissipative, 
partial differential equations (PDEs) are frequently used 
to study natural phenomena in many areas of physics. 
Generic PDEs such as the Kuramoto-Sivashinsky equa- 
tion and the complex Ginzberg-Landau equation have 
trajectories that are first exponentially attracted to a 
finite-dimensional inertial manifold and then settle in to 
a smooth global attractor. It is believed that a study of 
the covariant vectors or modes may lead to the identifi- 
cation of the inertial manifold [31^. The role played by 
the step structure in the Lyapunov spectrum for PDEs 
appears to be associated with convergence to the iner- 
tial manifold whereas for particle systems it is associated 
with the mechanically conserved quantities of the dynam- 
ics. Despite this difference in phase variables the angle 
distributions between covariant Lyapunov modes in both 
particle systems and PDEs seem to show some evidence 
of universal behaviour. 

The study of the covariant (C) Lyapunov vectors is 
less well understood but much progress has been made 
[271 [29H32] . Covariant Lyapunov vectors [28 follow the 
physically important features of the nonlinear dynam- 
ics as they align with the instantaneous directions of 
the stable and unstable manifolds. Thus instantaneous 
covariant vectors give information about the angles be- 
tween various manifolds in the multidimensional phase 
space dynamics, either angles between two different sta- 
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ble manifolds, two different unstable manifolds, or be- 
tween stable and unstable manifolds. 

In previous studies of the GS Lyapunov vectors an en- 
tropy based localization order parameter [22 was used as 
a probe of the structure of the vectors, and some aspects 
of the structure of the Gram- Schmidt vectors were ob- 
tained. Here we look more closely at the localization of 
the covariant vectors on average, the probability distri- 
bution of the localization and particular configurations 
which lead to features in both the distributions and in- 
stantaneous values of the order parameter. 

The manuscript is organized as follows. An introduc- 
tion to the method of calculating the covariant Lyapunov 
vectors. The description of the QOD model and the GS 
results for the Lyapunov spectrum and the Lyapunov 
modes. Some mention of the scalar properties of the 
system. The bulk of the paper introduces the Localiza- 
tion order parameter and discusses in detail the results 
which follow; the average localization, distributions of 
localization, configurations which lead to strong localiza- 
tion, instantaneous localization, strong localization and 
the scaling of localization. Next we look in detail at the 
angle distributions between all covariant Lyapunov vec- 
tors in the system and discuss the possibilities of tan- 
gencies. We finish with a review of our semi-analytic 
knowledge of the behaviour of the covariant Lyapunov 
modes based on an understanding of the dynamics of the 
GS Lyapunov modes. 

II. COVARIANT LYAPUNOV VECTORS 

The Benettin scheme for calculating the Lyapunov 
spectrum consists of time evolving a set of orthogonal 

(i) 

basis vectors which span the tangent space. Here j 
is the number of the basis vector and n — 1 is the step 
number (or discrete time). These vectors are periodi- 
cally re-orthogonalized using the Gram- Schmidt proce- 
dure and the exponential growth and decay rates of each 
vector give the Lyapunov exponents. This set of orthog- 
onal vectors are referred to as the Gram-Schmidt (GS) 
Lyapunov vectors. 

A number of methods for calculating covariant Lya- 
punov vectors have recently been proposed [24, 26 but 
the method we follow here is that due to Ginelli et. al. 
[25^, see also [29 for more details. Making the normal- 
ized Gram-Schmidt vector g^lu at phase point x^_i, the 
j^^ column of the matrix Gn-i^ the forward time evolu- 
tion in tangent space involves two steps. First, evolve the 
system forward in time using the tangent space dynamics 
Jn-i SO that = Jn-ig^-n ^ud sccoud, orthonormal- 
ize using a QR decomposition. So in matrix form the 
push forward in time is 

Gn = Jn-lGn-l (1) 

and the QR decomposition is 

Gn = GnRn- (2) 



where the columns of Gn are orthogonal. 

The j^^ covariant vector Vm is contained in the sub- 
space spanned by the first j GS vectors, so the matrix 
Vm which has covariant vectors as its columns, can be 
written as 

= GniCm (3) 

where Cm is an upper diagonal matrix of coefficients. 
The time evolution of Vm is constructed from the time 
evolution of Cm using Eqns ([T]), ([2| and (|3| so 

Vm — Gm-Rm Cm — Jm—lGm—lRm Cm 

= Jm-iym-lC^^R^^Cm (4) 

The key step is that if we choose C^R'^Cm — 
then time evolution of V becomes covariant as Vm = 
Jm-iVm-i and C~^R^^Cm = I implies that the time 
evolution of the C matrix is generated backwards in time 
by the R~^ matrices using 

Cm-l = RniCm (5) 

It is usual to normalize the columns of the coefficient 
matrix C at each step so that the covariant vectors are 
unit vectors. 



III. THE QUASI-ONE-DIMENSIONAL SYSTEM 

The model we consider consists of N hard disks (typi- 
cally 80) in a two-dimensional rectangular space L^x Ly. 
When Ly < 2cr, where a is the diameter of the disks, the 
space becomes a Quasi-One-Dimensional (QOD) system 
|5] and the hard disks remain ordered in the x-direction. 
For the system considered here we use Ly = 1.15cr, and 
scale Lx according to the density p = Ncj'^ / {L^Ly) with 
hard-wall boundary conditions in the x-direction and 
periodic boundary conditions in the ^/-direction ((H,P) 
boundary conditions) see Fig. [l] The phase space vector 
r = {q^p) consists of all particle positions q = (qi, qat) 
and momenta p = (pi,...,p7v) so the time evolution is 
composed of two components; the free-flight of particles 
between collisions and the change in momentum at col- 
lision 0. The significant advantage of using the QOD 
system is that both the Lyapunov exponents and the 
Lyapunov modes of the system can be obtained to high 
accuracy by standard numerical schemes with fast con- 
vergence rates. The exponents and mode structure for 
this system is well known consisting of a zero sub-space 
(where A = 0), and a number of transverse (T) and 
longitudinal-momentum proportional (LP) modes asso- 
ciated with the smallest positive and negative exponents 

The particle dynamics is symplectic [17] so the Gram- 
Schmidt (GS) Lyapunov vectors 5V = {Sq,6p)^ exhibit 
a characteristic structure. If the vector corresponding to 
the positive exponent is given by 
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FIG. 1. A schematic diagram of the Quasi One Dimensional 
(H, P) system with hard wall boundaries in the x-direction 
and periodic boundaries in the ^/-direction. The shaded disks 
are within the QOD system and the unshaded disks are the 
first periodic images above and below each particle. 



then the Lyapunov vector for the corresponding conju- 
gate negative exponent is given by 

^r4^+i-i = ("4^)- 

These conjugate relations are obeyed by the GS vectors 
and as they imply orthogonality, they do not apply to the 
covariant vectors which are in general not orthogonal but 
follow the directions of the stable and unstable manifolds. 
In the mode space the GS vectors occur in conjugate 
pairs and the majority of the dynamics occurs in the 
sub-space spanned by the pair of conjugate vectors 
This together with a detailed knowledge of the tangent 
space dynamics allows us to quantitatively understand 
the angle distributions between conjugate modes [9^. 

IV. LYAPUNOV EXPONENTS AND MODES 

It is well known O [TT] that the Lyapunov spectrum for 
the QOD system has exponents that appear in steps and 
associated with each of these exponents is a Lyapunov 
vector (or mode) that is delocalised and has one of two 
possible structures, either a transverse (T) mode or a 
longitudinal-momentum proportional (LP) mode. Occa- 
sionally, when the exponents corresponding to a T mode 
and a LP are equal (or approximately equal) the modes 
mix and become a time dependent linear combination of 
both mode forms (we refer to this as mode mixing) [6 . 
However, two features of the spectrum remained unex- 
plained. First the size of the step region in the expo- 
nent spectrum, and second the order in which the modes 
appear for a particular density. Recent calculations for 
systems of 100, 150 and 200 particles suggest that the 
thermodynamic limit for the spectrum has been reached 
at = 100 and that a fixed fraction of the spectrum con- 
sists of steps (and thus modes) [8 . The order in which 
modes appear can also be determined from Fig. [3) For 
any fixed density, the exponents when plotted as a func- 
tion of T mode number or LP mode number are to a good 
approximation linear with slopes that depend upon the 
density, as shown in Fig. [2] for a density of p = 0.5. Once 
the slopes of these lines are determined for a particular 
density by reading them off Fig. [3| graphs like that in 
Fig. [2] can be constructed and the T and LP exponents 



read off. The order of the steps is then determined by 
simply ordering the exponents from smallest to largest. 

0.25 




FIG. 2. (color online) The values of the Lyapunov exponents 
for the LP (blue) and T (red) modes of a system of = 
80 hard disks at a temperature of T = 1 and density p = 
0.5 plotted as a function of mode number. To a very good 
approximation the results are linear in mode number for both 
types of modes, with in this case a slope of 0.037654 for T 
modes and a slope of 0.031295 for LP modes. 

The slopes of the linear fits for the T and LP modes 
are shown in Fig. |3]as a function of density for 0.0003 < 
p < 0.8. We see immediately that for densities below 0.3 
the order of the steps is fixed, first an LP mode then a T 
mode, and then that pattern repeats throughout the step 
region. In fact this region of densities p < 0.1 has already 
been termed the kinetic region due to the behaviour of 
the largest exponent and the localization [22 . 

The functional form of GS and covariant modes are 
largely the same for the modes corresponding to posi- 
tive exponents n > but the modes corresponding to 
negative exponents are different. The negative covari- 
ant vectors are not orthogonal to their conjugate vectors 
and therefore there is a well defined linear combination 
of modes which they comprise. 

A. The Lyapunov Exponent Spectrum 

A systematic approach to the study of the Lyapunov 
spectrum as a function of density begins by looking for 
scaling behaviour in the exponents. For the QOD system 
it is known that the largest exponent at densities below 
p ~ 0.1 depends on the density as apln{/3p) as shown 
by Krylov [181 and this has been confirmed in previous 
simulation studies for the QOD system [22j and more 
generally for two-dimensional systems of disks f36]. As 
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FIG. 3. (color online) The slopes of the linear fits in Fig. [2]for 
systems of = 80 hard disks at a temperature of T = 1, as 
a function of density p. For the low density or kinetic region 
the slope for the LP (blue) and T (red) modes are almost the 
same over a large density range suggesting stronger scaling 
properties. However, at p = 0.3 the slopes cross over so that 
the LP mode has the largest slope. 



the Lyapunov spectrum for this system has been studied 
before we limit our description to features not previously 
observed. In Fig. Q, by plotting A/p against exponent 
number, we observe that for a large proportion of the 
Lyapunov spectrum in the kinetic region^ that is, expo- 
nent numbers between 40 to 160, Lyapunov exponents 
scale with density A ex p for the A/" = 80 disk system. 
This is observed in the whole positive half of the spec- 
trum, except for the first 40 exponents (for = 80), 
and includes all of the step region where the correspond- 
ing Lyapunov vectors are modes. It is known that the 
largest exponent scales as — p In p in the kinetic region so 
it is reasonable to assume that all exponents that are not 
in the linear region (that is 2 to 40) will scale in the same 
way as the largest. The proportion of the spectrum com- 
posed of steps and Lyapunov modes appears to remain 
fixed as the system size increases [8 so we expect that 
the scaling of the KS entropy is dominated by the scaling 
behaviour of the first 40 exponents. 



V. SCALAR PROPERTIES 

The thermodynamic and dynamical behaviour of hard 
disks in narrow channels has been studied before [38] as 
a function of Ly. This system becomes QOD when Ly < 
2 and it is found that the collision frequency and the 
Lyapunov spectrum scale strongly with the width of the 
system Ly. This implies a strong scaling of both the 
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FIG. 4. (color online) The scaled Lyapunov spectrum (A/p) 
for the positive exponents in the kinetic regime for systems 
of A/" = 80 hard disks at a temperature of T = 1 (that is, 
densities ranging from p = 0.0003 to p = 0.03). The modes 
occur in the same order throughout this regime, with an LP 
mode followed by a T mode. Notice that for all densities 
below p = 0.1 the values of \n/p are the same for exponent 
numbers 40 < n < 160. In contrast, the scaling is not linear 
in density for the first 40 exponents. 



largest exponent and the KS entropy with Ly. 

If a system of A^ particles has A^^ collisions in time 
t then the single particle collision frequency is given by 
V2 = 2nt/{Nt). For the QOD system the single parti- 
cle collision frequency is directly related to the potential 
contribution to the pressure as 



V2 



2{Pv - kT) 



(8) 



where Ap is the average momentum transfer per collision. 
Here we see that the single particle collision frequency at 
a density of 0.01 changes by a factor of 3 as Ly goes from 
1 to 2, that is over the range of Ly values that correspond 
to the QOD system. 

There have been a number of theoretical studies of the 
KS entropy per particle of hard-particle systems in two 
and three dimensions [3T. In the low density limit the 
KS entropy is expected to be linear in the single-particle 
collision frequency V2 given by 



h 



KS 



N 



V2A{-\np^B^0{p)) 



(9) 



where A and B are constants, and will thus change dra- 
matically with the width of the QOD system Ly. As 
the KS entropy is simply the sum of the positive Lya- 
punov exponents its scaling behaviour is determined by 
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the dominant scaling of individual groups of exponents, 
in particular the first 40 exponents for a 80 disk system, 
As we have only considered one fixed value of Ly — 1.15, 
we have not attempted a discussion of the strong scaling 
of these properties with Ly . 



VI. LOCALIZATION 

A. Order Parameter 

The contribution of particle i to the j^^ Lyapunov vec- 
tor Sr^^^ can be written as x?^ = + (^P?^) • 
By definition is positive and as the Lyapunov vec- 
tor is normalised, Xi"^^ ^ 1 foi" each particle. Therefore, 
for each Lyapunov vector j the set {Xi^^} are positive 
random variables which sum to one. It is natural to con- 
struct an entropy function for each Lyapunov vector [22] 

from the set {Xi"^^} and this gives a measure of the local- 
ization of the Lyapunov vector 

W^i/)W = lexp(^-f^X?^lnx?^). (10) 

This measure W^^^ {t) is close to zero for a highly localized 
vector and close to one for a strongly delocalized vector. 
It gives an instantaneous measure of the localization for 
each vector j and represents the proportion of the N 

particles that contribute to the Lyapunov vector. Indeed 

(i) 

NWj^ could be considered the participation number - the 
number of contributing particles - which is independent 
of N if the vector is strongly localized on a subset of 
the particles. In this way it relates more directly to the 
inverse participation ratio I2 |19H21) which is also used as 
a measure of the number of effective degrees of freedom. 

(i) 

However, no simple exact relation between NWj^ and 
Y2 can be given. 

In what follows, we will investigate all aspects of the lo- 
calization for the QOD system - the average localization 
of both GS and covariant vectors, the probability dis- 
tribution for the localization, configurations which lead 
to strong localization, instantaneous localization and the 
connections between localization and angle distributions. 



B. Average Localization 

The average localization for a QOD system of = 80 
disks at two different densities is given in Fig ([5|. One 
density p = 0.003 is typical of the kinetic region and the 
other p = 0.8 is typical of the dense liquid phase. It has 
been observed [22 that the localization of the GS vectors 
as a function of density generally takes one of two forms: 
one typical of a low density kinetic region, and the other 
typical of a dense liquid. Strong localization occurs in the 



kinetic region when the density is less than 0.1, and here 
a group of Lyapunov vectors associated with the largest 
exponents show small values of wj^\ Indeed, at very 
low densities a number of vectors all show a similar form 
of strong localization which appears to depend linearly 
on the vector number [23] in the localization spectrum. 
For the GS vectors the conjugate relation Eqn ^ implies 
that the localizations for the two conjugate vectors are 
identical and this is indeed what is observed in Fig. ([5|. 
The behaviour of the localization for the range of vectors 
140 < j < 160 (depending on system size) is associ- 
ated with the different localization values for T modes 
and LP modes. The GS localization spectrum at high 
density {p = 0.8) is essentially exponential in the vector 
number and more delocalized than the covariant vectors 
which look linear in the vector number. At low density 
{p = 0.003), as a function of vector number, the local- 
ization is initially linear before returning to the usual 
exponential shape. In Fig (|5| we see the same effect at a 
density of p = 0.003 where both the GS and covariant lo- 
calizations depend linearly on the vector number between 
1 and about 35. The randomly dropped brick model of 
[23] gives an estimate of the number of most localized GS 
vectors that can co-exist and not violate the orthogonal- 
ity condition. For 80 disks it is possible to have 40 such 
vectors but on average random positioning reduces this 
to approximately 35. This is in good agreement with the 
results reported in Fig. ([5|. 

At high density the covariant modes with positive ex- 
ponents have similar localizations to those associated 
with negative exponents, but at low density there are 
small differences. The localization of the Lyapunov vec- 
tors associated with the most negative exponents appears 
to be different to the localizations of the conjugate vec- 
tors, suggesting the possibility of some convergence is- 
sues. The Gram-Schmidt vectors, from which the covari- 
ant vectors are calculated, satisfy an exact conjugate rela- 
tion for the localization however, the linear combinations 
which make up the covariant vectors need not. However, 
the time reversal symmetry of this system suggests that 
the localizations for positive and negative vectors should 
be the same. 

The covariant vectors at low density have a localization 
spectrum that is again linear but with smaller slope (al- 
most constant). For covariant Lyapunov vectors there is 
no longer a strict orthogonality constraint so subsequent 
Lyapunov vectors can contain a significant component of 
the first, second or higher, Lyapunov vectors. Therefore 
the localization characteristics of the first covariant vec- 
tor can infiuence the localization of subsequent covariant 
Lyapunov vectors more strongly. The j^^ covariant vec- 
tor can be written as a linear combination of the first j 
GS vectors as Vj = X]i=i where q are coefficients, so 
for Vj to be localized we require all of ^1, to be lo- 
calized on essentially the same group of particles. Even if 
all component vectors are strongly localized it does not 
necessarily follow the the covariant vector is localized, 
although the numerical evidence suggests that the same 
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FIG. 5. (color online) The average localizations W for both 
the GS (red) and the covariant (blue) vectors of a QOD sys- 
tem of 80 hard disks at a temperature of T = 1 , at two differ- 
ent densities. In the first figure the density is p — 0.003 and 
for the second figure the density is p = 0.8. There is clear evi- 
dence of strong localization in the low density state for vector 
numbers j < 35. In the mode region the localization splits 
into two branches, one for the T modes and the other for the 
LP modes. 



number of strongly localized vectors occurs in both the 
GS basis and the covariant vectors, the covariant ones 
are on average more localized than the GS vectors. 



C. Localization Distributions 

The localization W^\t) is defined at each instant in 
time so we can construct a probability distribution for 
this quantity for each Lyapunov vector in Figs. ([6| and 
0. This gives more detailed information than the av- 



erage localization < W^'' > in Fig. ([s]). The average 
localization as a function of density shows that the zero 
modes and the T and LP modes are only weakly density 
dependent. The major feature of the localization in the 
continuous part of the spectrum is that the vector cor- 
responding to the largest exponent becomes more and 
more localized with decreasing density or increasing N. 




FIG. 6. (color online) The probability distribution of the 
localization W for the covariant vectors of a QOD system of 
80 hard disks in the kinetic regime at a density p — 0.003 
and temperature of T = 1. The axis labeled 0.0 to 1.0 is 
the value of the localization, the axis labeled to 127 is the 
vector number and the vertical axis is the probability. The 
most probable localization value for the first 35 vectors is 
the minimum localization of 2/N. There is also a small cusp 
at 4:/N corresponding to a linear combination of two most 
strongly localized configurations. For vectors larger than 40 
the localization distributions are smooth and very similar to 
those of higher density systems where the distribution has a 
smooth maximum. 



In Fig. (|6| the distribution of localization for the first 
127 Lyapunov vectors at a low density {p = 0.003) for 
a QOD system of 80 disks are shown. Initially there is 
a dominant contribution at the position of the most lo- 
calized configuration W = 2/N with a barely visible blip 
at W = A/N corresponding to four equal contributions 
to the set {Xi}- However, with increasing vector number 
the most localized peak diminishes apparently linearly 
and the distribution becomes broad and smooth. For 
Gram-Schmidt vectors this most localized contribution 
occurs only for the first quarter of the expanding Lya- 
punov vectors, but for the covariant Lyapunov vectors 
the influence of this most localized contribution can con- 
tribute to more vectors as seen in Fig. (|6|. Therefore 
the structure of the most localized configuration plays 
an even more important role for the covariant Lyapunov 
vectors. At high density there is a remnant of the strong 
localization configuration evident in Fig. ^ but its am- 
plitude is only 1/3 of that at low density and decays more 
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FIG. 7. (color online) The probability distribution of the 
localization W for the covariant vectors of a QOD system of 
80 hard disks at a density of p = 0.8 and temperature of 
T = 1. The figure is to be compared with the kinetic region 
in Fig ([6| but here the vertical scale is one third smaller. 
Otherwise the axis labels are the same as Fig ([6|. There IS a 
small remnant of the most localized configuration in the first 
few vectors. 



quickly than linearly with vector number, and with the 
scale change in the graph the distribution with a smooth 
maximum is more apparent. 



D. Strongly Localized Configurations 

In Fig. ([6|, for each of the first 35 Lyapunov vectors, 
there is clear evidence of an instantaneous localization of 
W = 2/N = 0.025 which corresponds to a configuration 
with X = 0-5 oil two neighbouring particles remaining 
for long periods of time. The position of the localized 
pair of particles is not important but in the transition of 
the localization from one pair of particles to another a 
variety of things can happen. The first possibility is that 
the localization on two particles decays from an initial 
value of 1/2 to zero as the localization on a different two 
particles increases from zero to 1/2. This behaviour is 
clear in the first GS vector (in Fig. ([6|) between collision 
numbers 80 and 90 where there is a plateau value of 4/A^ 
as predicted above. There are other possibilities as we 
see localization values above A/N for short periods which 
involve more than 4 particles, and other smaller blips as 
the pair of vectors moves one place in either direction. It 
is evident that at this time the second covariant vector is 
essentially equal to the first GS vector (although it must 
contain some small component of the second GS vector). 

The most localized configuration is when = 0-5 
for both i and i -\- 1 and is zero for all the other values 
of i (see Fig. Is] (a)), then the localization measure is 




given by 



2/N. Clearly the vector is localized 



FIG. 8. (color online) (a) The configuration which gives the 
strongest localization W^^"^ for the QOD system of hard disks. 
Here Xi'^'^ =0.5 for i = 17 and i = 18 and thus this configu- 
ration gives W^-^^ = 2/N. (b) A linear combination of two of 
the most localized configurations gives a range of values of lo- 
calization depending upon the particular linear combination. 
If Xi^'^ = 1/a for i and i -\- 1 and Xk"^ = 1//3 for /c and k -\- 1 
the localization depends on the values of a and /3 which are 
not independent. When a — 0, wj^-* = 2/N and when a — /3 
then W^^^ = 4/N. 



on only two particles so that W^^"^ follows a strict power 
law in the number of particles A/", with the power equal 
to —1. Another possible strongly localized configuration 
is a mixture of two of these most localized vectors (see 
Figjs] (b)). As this is a mixture, we expect that Xi — 

1/a for i and z + 1 and for x^k^ — 1//^ ^ k -\- 1 
where P = a/{a — l). The localization for this mixture is 
]^0)(t) = a(a—l)^"'^/"~"'^Y^ which has a maximum when 
a = 4. In fact this combination only gives localization 
values between 2/N and A/N. 



E. Instantaneous Localization 

In Fig. ^ for N = 80, the low density average lo- 
calization for the covariant Lyapunov vectors shows an 
initial strongly localized region associated with the first 
35 Lyapunov vectors. A similar behaviour has been seen 
previously in the average localization of the GS vectors 
and there it was argued that each of these vectors was 
dominated by the most localized configuration in Fig. (|8| 
and the required orthogonality limited the number of vec- 
tors. Here we know that the covariant vectors are linear 
combinations of the GS vectors, but that in general the 
covariant vectors are more localized than the equivalent 
GS vectors. Here we look at the dynamics of the localiza- 
tion in order to understand this property of the covariant 
vectors. In Fig. ^ we show the evolution of the for 
each particle, from the GS vectors, over 100 collisions 
and at the same time the evolution of the localization or- 
der parameter W{t). For vector 10 there are two places 
where we see the most localized configuration around 60 
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and 90 collisions but there are larger excursions that do 
not correspond to either of the two strongly localized con- 
figurations. Similarly, for vectors 20 and 30 there seem 
to be no most localized configurations and increasingly 
larger excursions in WN{t). 
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FIG. 9. (color online) The localization contributions from 
each particle Xi'^^ calculated as a function of collision num- 
ber for a selection of the GS vectors of a QOD system of 
= 80 hard disks at low density p = 0.003 and a temper- 
ature of T = 1. The upper panels are the values of x?'* 
for each particle for Lyapunov vectors j — 10, 20, 30 with the 
lower panels giving the instantaneous value of the localization 
for the system W{t). The vertical axis on the upper panels 
is the particle number and on the lower panels is the local- 
ization. The horizontal axis is the collision number for all 
graphs. The color changes on the upper panels correspond 
to levels 0.05,0.15,0.25,0.35 and 0.45 and the corresponding 
colors are in order: green, white, blue yellow, orange and red. 
The vector j — 10 is strongly localized and for several periods 
its value is 2/N = 0.025, so its configuration corresponds to 
Figure (|8|. The localization is less strong for vectors j — 20 
and j = 30. 



In Fig. (10) we show the same time evolution of Xi^ 
for each particle, for the covariant vectors, over the same 
100 collisions. Vector 10 is initially localized on two par- 
ticles for the first - 10 colhsions so W{t) = 2/N = 0.025 
and then quickly splits into two two particle pairs where 
W{t) ~ 0.05. This remains until ~ 75 collisions when the 
localization returns to the value W{t) = 0.025. The other 
two vectors 20 and 30 also show short regions where the 
strong localization W{t) = 0.025 occurs but for vector 
20, W{t) ~ 0.05 is more common and then at vector 30 
larger excursions from the two most localized configura- 
tions involving more than 4 particles are observed. 

The time evolution of x''/^ foi" each particle, for the co- 
variant vectors at a high density of p = 0.8, over 100 col- 
lisions shows that none of the vectors show the strongest 
localization but there are periods where the localization 
is principally on four particles. 




FIG. 10. (color online) The localization contributions from 
each particle Xi^'^ calculated as a function of collision number 
for a selection of the covariant Lyapunov vectors of a QOD 
system of = 80 hard disks at low density p = 0.003 and 
a temperature of T = 1. The upper panels are the values 
of Xi'^'* for each particle for covariant Lyapunov vectors j — 
10, 20, 30 with the lower panels giving the instantaneous value 
of the localization for the system W{t). The axes and color 
changes are same as in Fig. (|9|. There is evidence of strong 
localization in both j = 10 and j — 20, but for j — 30 the 
localization is less strong. 



F. Strong Localization 

To establish asymptotic strong localization in the ther- 
modynamic limit we need to look at the localization as a 
function of N . If the system exhibits asymptotic strong 
localization then only a finite number of degrees of free- 
dom are involved and this does not change with system 
size N . Localization on only two particles is the most 
localized configuration possible. Another common situa- 
tion is when the localization is confined to a finite number 
of particles, strictly less then N . In that case the entropy 
in the exponent of Eqn (10) is a constant so the partic- 



ipation number NW^^ (t) which is constant independent 



of N. In Fig. (11) we plot the participation number for 
QOD systems of N = 20, 40, 80, 160 particles, and the re- 
sults show that the region of strong localization already 
identified in the average localization figures, is indeed 
strongly localized. 



G. Localization Scaling 

Earlier we saw that in the kinetic regime the Lyapunov 
spectrum exhibited a scaling relation with the exponents 
being linear in density. Here we see that in the same 
range of exponent numbers, that is 45 to 160 for a 80 
disk QOD system, the average localization is a constant 
independent of density. The results are shown in Fig. 
(12), and the inset in that graph shows how the initial 



45 exponents do change with density. 



9 



NW(N) 




0.8 



w 




FIG. 11. (color online) The participation number N (w^-*^ 

for the covariant vectors of a QOD system of hard disks 
as a function of the number of particles N at low density 
(p = 0.003) and a temperature of T = 1. The horizontal axis 
is a continuous variable constructed from the vector number j 
using X = {j — l/2)/2N. There is evidence of strong localiza- 
tion in the first 20% of the covariant vectors corresponding to 
positive exponents. As the system size is increased the same 
proportion of vectors remains strongly localized. 



VII. ANGLE DISTRIBUTIONS 



FIG. 12. (color online) The localization spectrum (W) for the 
positive covariant Lyapunov vectors for systems of 80 hard 
disks at a temperature of T = 1 with densities in the kinetic 
regime ranging from p = 0.0003 to p = 0.03. The curve 
for 0.0003 is red, 0.001 is blue, 0.003 is green, 0.01 is black 
and 0.03 is brown. Above 140 the modes separate into two 
branches and alternate between the values for LP modes then 
T modes for all densities, and again the localization of each 
mode is independent of the density. For densities below 0.1 
the values of W are the same for vectors in the range 45 < 
n < 160. The first 50 Lyapunov vectors show some small 
systematic variation in localization as evidenced in the inset. 
Here the localization becomes smaller with smaller density. 



The distribution of angles between covariant Lyapunov 
vectors j and j -\-dis observed to be largely independent 
of j but strongly dependent on d. A simple geometric 
argument has been proposed to explain this observation 
[39 . Imagine placing an initial unit vector randomly di- 
rected at the origin of a high-dimensional space. Placing 
a second vector at random can be achieved by choosing a 
circle of unit radius containing the head of the first vec- 
tor, and choosing a random point on this circle as the 
end-point of the second vector. The angle distribution 
between these two vectors is uniform. For the third vec- 
tor we consider a sphere containing the heads of vectors 
1 and 2, and choose a random point on the sphere to be 
the head of the third vector. The resulting distribution 
is p{Ois) = sin(^i3). It follows that the distribution of 
angles between two vectors with indices differing by d is 
p{Od) = s'm^^~^\Od)' Here we see that this geometric ar- 
gument works reasonably well at high densities, but at 
low densities in the kinetic regime, the angle distribu- 
tions are affected by the strong localization of some of 
the covariant Lyapunov vectors. Indeed, looking back at 
the high density results there still remains a small signa- 
ture due to localization in the first few Lyapunov vectors. 
The angle Oij between a pair of Lyapunov vectors i and j 



is calculated at each collision and a distribution function 
for the angle constructed. 



In Fig. (13) we show the angle distributions between 



nearest neighbour vectors for a low density state in the 
kinetic regime for all pairs of vectors. The probability 
distribution for the angle is at first glance uniform 

but with noticeable peaks near parallel (0, tt) and orthog- 
onal (7r/2). This system shows strong localization in the 
first 50 or so Lyapunov vectors and the dominant first 
vector makes a major contribution to all the subsequent 
strongly localized vectors. This effect leads to the peaks 
for near parallel configurations of nearest neighbour pairs 
of vectors for i = 1, .., , 50, becomes smaller with increas- 
ing i and is absent after about z = 50. The peaks near 
and TT suggests that pairs of vectors come arbitrarily 
close, almost parallel, in both the unstable manifold and 
the stable manifold. There is a similar central orthogonal 
peak at 7r/2 which decays very quickly with increasing i 
and disappears completely after z = 25 or 30. While it is 
clear that a generic chaotic trajectory cannot develop an 
exact tangency, there appears to be a possible intermit- 
tency mechanism that inflates the measures near and 
TT. In the centre of Fig. (13) between vectors i = 160 to 
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i = 200 which correspond to the modes, we see a strong 
preference for close to orthogonal configurations but also 
a lone peak. Apart from this region, Fig. (13) looks 
remarkably symmetrical with the distribution of 
almost to equal to the distribution for conjugate vector 
pair ^4Ar-i,4Ar-i+i. 




3 "'0 



FIG. 13. (color online) The probability distribution for the 
angle between nearest neighbour covariant Lyapunov 

vectors for a system of 80 QOD hard disks at a density p = 
0.003 and temperature of T = 1. The axis labeled to 3 is 
the angle in radians, the axis labeled to 300 is the vector 
number and the vertical axis is the probability. The angle 
distributions initially show a strong preference for angles of 
and TV and a smaller preference for 7r/2. These preferences 
decay slowly with increasing vector number and disappear 
completely after vector number 50. 



In Fig. (14) we show the angle distributions between 



second nearest neighbour vectors for a low density state 
in the kinetic regime for all of the possible second neigh- 
bour pairs. Here we see similar effects in the distribution 
due to strong localization of the covariant Lyapunov vec- 
tors. The peaks in the first 50 pairs of vectors at and tt 
are smaller than in Fig (13) while the peak at 7r/2 is a lit- 



tle larger and remains apparent for longer. There is again 
evidence of strong structure in the mode region which is 
largely orthogonality between second neighbour vectors 
but there are two peaks with positions much closer to 
and TT. Again there is symmetry between conjugate 
pairs of vectors angles with the distributions for 
and 04N-i,4N-i-\-i the same. 

As we move from 2nd neighbours to 3rd neighbours, 
4th neighbours and beyond, the central peak at 7r/2 in- 
creases and the peaks at and tt decrease. 

In Fig. ( 15 ) we give a typical higher density result for 
the angle distribution of all 319 nearest neighbour pairs of 
Lyapunov vectors at p = 0.8. While the angle dis- 

tributions are approximately constant in the continuous 
part of the spectrum there are features associated with 
remnants of the low density localization of individual vec- 
tors present at the beginning and end of the spectrum. 




FIG. 14. (color online) The probability distribution for the 
angle Oi,i-\-2 between second nearest neighbour covariant Lya- 
punov vectors for a system of 80 QOD hard disks at a density 
p = 0.003 and temperature of T = 1. The axes labels are 
the same as Fig. (13). Here the initial preference for and tt 



are less strong compared with the preference for 7r/2. Again 
these preferences decay and have disappeared completely af- 
ter vector number 60. 



and a preference for orthogonality for nearest neighbour 
vectors in the mode region. Otherwise the distributions 
show good symmetry for conjugate pairs of angles 
and ^47v-i,4Ar-i+i. In Fig. (16) we show the angle distri- 
butions between second nearest neighbour pairs of 
vectors in the high density regime for all pairs of second 
nearest neighbour vectors. Here the background distri- 
bution is approximately sinusoidal with deviations for a 
group of the first and last pairs of vectors, as well as 
strong structure in the mode region. The distributions 
show the same conjugate symmetry as before. 

Fig. (17) shows how the distribution for the first near- 
est neighbour angle 0i^2 varies with density for an 80 
QOD particle system. At the lowest density of 0.0003 



the peaks at 0, 7r/2 and tt dominate the distribution. 
Here the most localized configuration of a covariant vec- 
tor is the dominant effect, so two vectors which are of 
most localized type, but supported on different pairs of 
particles are clearly orthogonal and contribute at 7r/2. If 
they are supported on the same pair of particles then this 
will lead to parallel or anti-parallel situations which con- 
tribute to peaks at and tt. The parallel and anti-parallel 
configurations need to be more closely analysed as their 
nature cannot be assessed on the values of alone, but 
depend on the particle component contributions which 
lead to the values of Xp^- Also different contributions 
will occur when the two vectors share only one parti- 
cle in common. The distributions change systematically 
from low density, where the effects of localized covariant 
vectors dominates, to high density where the effects of 
localization have almost disappeared. 
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FIG. 15. (color online) The probability distribution f{Oi,i-\-i) 
for the angle between nearest neighbour covariant Lyapunov 
vectors for a system of 80 QOD hard disks at a density p = 0.8 
and temperature of T = 1. The axes labels are the same 
as Fig. (13). The three-dimensional plot of all vector pairs 



shows a high degree of symmetry between nearest neighbour 
angles corresponding to positive and negative pairs of vec- 
tors. Strong structure is also present in the angles between 
Lyapunov modes. Here we show results for pairs i = 1 to 319. 
At this density the angle distributions are only approximately 
constant, and the differences are associated localization. The 
distribution is slightly peaked at 0, 7r/2 and tt but this fea- 
ture disappears as the vector numbers increases from 1 or 
decreases from 319. 



VIII. TANGENCIES 

So far the angle distributions we have considered are 
the angles between close neighbouring covariant vectors 
that are almost always contained in the same manifold, 
either both in the unstable manifold or both in the sta- 
ble manifold. To consider the possibility of tangencies, 
which we interpret as the finite probability of an angle 
of or TT between one vector contained in the unstable 
manifold and one vector contained in the stable manifold, 
we consider the distribution of angles between conjugate 
covariant vectors i and AN — This criteria has been 
proposed as a method to find the smallest angle between 
the stable and unstable spaces ^ and probe the break- 
down of hyperbolicity. 

In this QOD system we have seen that the smallest 
angles between stable and unstable covariant vectors oc- 
curs for conjugate pairs of vector modes, and typically 
the first conjugate pair of modes contains the smallest 
angle [9 . At high density the first mode is a conjugate 
pair of transverse (T) modes while at low density it is 
a conjugate pair of LP modes. The distribution of an- 
gles between all conjugate modes at high density is given 
in Fig. (18) and we see a clear distinction between the 



conjugate angle between modes and between vectors in 
the continuous region. There is a clear peak near zero 
for the first conjugate pair of T modes which is bounded 
away from zero and each subsequent pair of modes gives 



FIG. 16. (color online) The probability distribution for the 
angle ^i,i+2 between second nearest neighbour covariant Lya- 
punov vectors for a system of 80 QOD hard disks at a density 
p = 0.8 and temperature of T = 1. The axes labels are the 
same as Fig. ( |13[ ). At this density the distributions are very 
similar except for the first 10 or so vectors. This is very simi- 
lar to the systematic differences noted in Fig (15) in the first 
few vectors. 



f(e)de 



0.01 



0.001 




FIG. 17. (color online) The natural logarithm of the proba- 
bility distribution for the angle ^1,2 for the first nearest neigh- 
bour pair of Lyapunov vectors for a system of 80 hard disks at 
temperature T = 1 as a function of density. At very low den- 
sity, in the kinetic region, the distribution is sharply peaked 
at 0, 7v/2 and tt but as the density increases the peaks be- 
come less sharp and smaller, until at p = 0.8 the distribution 
is almost constant in angle. Using a logarithmic scale on the 
vertical axis accentuates the differences between distributions 
at the lowest densities. 
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FIG. 18. (color online) The probability distribution for the 
angle ^i,4Ar-i+i between conjugate pairs of vectors at a density 
of p = 0.8 and temperature T = 1 for a system of = 80 
QOD hard disks. Here the axis labelled 1 to 160 is the vector 
number for the conjugate pair of modes. These begin at 1 
and 2 with the zero modes and increase so that 160 is the 
angle between the vectors corresponding to the largest and 
smallest exponents. In the mode region there are strong peaks 
but these are bounded away from and tt. The continuous 
spectrum vectors all show the same orthogonality. 




FIG. 19. (color online) The probability distribution for the 
angle Oi^AN-i+i between conjugate pairs of vectors at a density 
of p = 0.003 and temperature T = 1 for a system of iV = 80 
hard disks. This shows similar systematic changes with mode 
number compared with the high density result. There is a 
feature in the angle distribution associated with the strongly 
localized vectors in the positive and negative halves of the 
spectrum. 



p=0.003 



a different peak, all of which are closer to it/ 2. The vec- 
tors that are not modes all show very good orthogonality, 
strongly peaked at 7r/2. 

At lower density, p = 0.003, see Fig. (19), the distri- 



bution between the first conjugate pair of modes (here 
LP modes) approaches and tt but the distribution goes 
to zero at both ends, so again no possibility of tangency 
although the distribution is no longer strongly bounded 
away from tt. Indeed the distributions for all pairs of 
are much broader and less peaked than at high density. 
However, here there is a clear gap and no tangencies are 
possible. The vectors in the continuous part of the spec- 
trum all show strong orthogonality with dominant peaks 
at 7r/2. 

The angle distribution between the first conjugate LP 
modes at lower density {p = 0.003) are shown in Fig. 
(20). Although these distributions are the 3^^ line m 
Fig. (19), we have enlarged the vertical scale so that 



it is clear that the distributions are bounded away from 
possible tangencies at and tt, despite the fact that the 
distributions are broad and largest close to these limits. 

Previous theoretical studies of the time evolution of a 
pair of conjugate modes [7 showed that the dynamics 
approximately separates from the dynamics of the whole 
tangent space, and can be developed for a pair of T modes 
by considering the product of 2 x 2 reduced R matrices 
that produce the backward time covariant evolution |9]. 
These matrices have the form 



Cm ^mCm 
Cm 



(11) 



where ("m is the local expansion rate for the m*^ cycle of 
free-fiight and collision, and Cm is related to the free flight 



f(e)de 




FIG. 20. (color online) The probability distribution for the 
angle between first two conjugate pairs of LP modes at a 
temperature 1 and density of p = 0.003, for a system of 80 
hard disks. Notice that the distribution is bounded away from 
(or tt) and goes to zero fast enough to avoid tangencies. 



time (usually a little smaller than Tm)- The cosine of the 
angle between these conjugate modes can be written as 



cos Oi 



f 



VP 



1 



(12) 



where the function / depends on the sets of values {Q} 
and {cj} contained in all the Rj matrices along the path. 
For an evolution of n — m steps backward from m the 
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function is 

m 

/(m,n,{Cj},{c^}) =Cn+l + Cn+i--C|-iCj. (13) 

In the limit as n ^ oo the vector converges to the covari- 
ant vector and as long as Q > 1, the angle Oij is bounded 
away from or tt. To see the boundedness explicitly as 
n ^ oo we set (j = ( and Cj = r for all j and then 
/ = r/(l — C~^) which is finite as long as 7^ 1- As long 
as / remains finite in Eqn. (12), then cosine of the angle 
is bounded away from 1 or —1. 



IX. CONCLUSIONS 

We have studied the Lyapunov spectrum and vectors, 
both GS and covariant, for QOD systems over the full 
range of physical densities. The first observation is that 
the exponents in the step region depend linearly on the 
mode number for both T and LP modes and given the 
slopes of these linear dependences as a function of den- 
sity, we can predict the order in which the modes appear, 
or equivalently the order in which the steps appear in the 
spectrum. Next we observe that both the Lyapunov spec- 
trum and the localization obey simple scaling relations in 
the kinetic regime for the same group of vector numbers 
- typically numbers 40 to 160 for an 80 particle system. 
Here the Lyapunov spectrum is linear in density and the 
localization is independent of density. 

The average localizations for covariant modes is gen- 
erally significantly smaller and thus more localized than 
the results for the GS modes. This is observed at low 
density in the kinetic region and also at high density. 

The strongly localized configurations for GS vectors 
consisting of a dominant contribution from two neigh- 
bouring particles are again apparent and the covariant 
vectors being linear combinations of GS vectors show the 
same feature. However, the removal of the orthogonality 



condition in going from GS to covariant vectors allows 
one strongly localized configuration to dominate a whole 
group of covariant vectors whereas it can only contribute 
to one GS vector. This is evident for the low density aver- 
age localization where many of the strongly localized vec- 
tors have essentially the same average localization, and 
is connected with the instantaneous localization. 

There is very clear evidence of asymptotic strong lo- 
calization in the kinetic regime for the first 20% vectors 
of the covariant spectrum where the participation ratio 
NW{N) is independent of system size. 

The simple scheme for determining the angle distribu- 
tions between neighbouring covariant vectors proposed 
here {p{Od) = sm^~^{Od)) is useful in understanding the 
underlying behaviour of the distributions but the results 
are affected by the degree of localization and this is deter- 
mined by the density. At low density, strong localization 
has a significant effect however, there is some strong lo- 
calization effect at all densities - limited to the first and 
last few vectors at high density. 

There is no evidence for tangencies in the QOD sys- 
tem. Indeed, if two vectors become parallel then the ma- 
trix Vm loses rank and the method will fail. This is never 
observed for this QOD system in finite time. The closest 
approach to tangency is for the vectors associated with 
the first positive and negative modes. Only at low densi- 
ties in the kinetic regime do the first vectors approach 
or TT. The angles between conjugate vectors in the unsta- 
ble manifold and the stable manifold show a high degree 
of orthogonality suggesting that two manifolds are gen- 
erally orthogonal. 
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